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This study considers the periodic orbital period of an n-body system from the perspective 
of dimension analysis. According to characteristics of the n-body system with point masses 
(m1, M2, ..., Mn), the gravitational field parameter, a ~ Gmim,;, the n-body system reduction mass 
Mn, and the area, An, of the periodic orbit are selected as the basic parameters, while the period, 
Tn, and the system energy, |E,,|, are expressed as the three basic parameters. By using Bucking- 
ham z-theorem of dimensional analysis, these two relations can be reduced to a dimensionless form, 
which can surprisingly produce only one dimensionless 7, respectively. Because there is only one 7, 
therefore the 7 must be a constant. Since the two-body system is a special case of the n-body, we can 
uniquely determine the two constants by using the two-body Kepler’s third law. Thus, the n-body 


system Kepler’s third law is deduced and is given by Te |e! 2 
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A numerical validation and comparison study was hence conducted. 
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INTRODUCTION 


One of the central and vivid problems of celestial me- 
chanics in the 18th and 19th centuries was about the 
description of the motion of the Sun-Earth-Moon sys- 
tem under the Newtonian gravitation field, depicted in 
Figure 1. Notable work was done by Euler (1760), La- 
grange (1776), Laplace (1799), Hamilton (1834), Liou- 
ville (1836), Jacobi (1843), and Poincaré (1889) [1] and 
Xia (1992) [2]. The study of the motion of two bodies 
was solved by Kepler (1609) and Newton (1687) early in 
the 17th century. For the elliptic periodic orbit of 2-body 
system in Fig 1, Kepler’s third law of the two-body sys- 
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the gravitation constant, G = 6.673 x 107 lim?kg7!s7?, 
the orbit period, T, the total energy of the 2-body sys- 
tem, |E|, and point masses mı and mg. 


where 


tem [3] is given by T|E|?/? = Gmime 


FIG. 1: The Sun-Earth-Mond system and 2-body system 


However, the 3-body system depicted in Figure 2, can- 
not be solved analytically because unlike the 2-body 
problem, the 18 variables that describe the system can- 


not be reduced to a single variable. Simplification of the 
two-body problem was allowed by invariance and con- 
served quantities as ‘first integrals’. It was proven that 
it is impossible to reduce the 18 variables of the 3-body 
problem in order to produce an analytic solution. 
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FIG. 2: 3-body system 


Notwithstanding that the analytic solution cannot be 
found, it is possible to find a numerical solution for the 3- 
body problem, in which the study of the periodic 3-body 
orbit has received particular attention in recent years [4— 
11]. The figure-eight orbit was discovered numerically in 
1993, using the principle of the least action by C. Moore 
[4]. Chenciner and Montgomery [5] rediscovered the so- 
lution by using the shape space least action principle in 
2000. In 2013 Suvakov and V. Dmitraginovié [6] made a 
breakthrough and produced epic results, they found 13 
new distinct collisionless periodic orbits of the Newtoni- 
an planar 3-body system with an equal mass and zero 
angular momentum. In 2017 Li and Liao [10], Li, Jing 
and Liao [11] reported their breakthrough new finding: 
695 periodic orbits of planar 3-body system with an equal 
mass and zero angular momentum, as well as 1223 peri- 
odic orbits of the planar 3-body problem with an unequal 
mass and zero angular momentum, respectively. 

Corresponding to the new finding of more and more 
periodic obits, a fundamental law was discovered by nu- 


merical experiments [9-11], which states that the 3-body 
system Mı, M2 M3) may obey a law, which is similar to 
the law of harmonies, name Kepler’s third law. Anal- 
ogous to the Kepler’s third law of the two-body prob- 
lem, Suvakov and V. Dmitrasinovié [9] proposed a gen- 
eralized 3-body Kepler’s third law: T|E|°/? = constant, 
where || denotes the total kinetic and potential energy 
of the 3-body system, where T is the period of periodic 
orbit. However, they pointed out that ”the constant on 
the right-hand-side of this equation is not “universal” in 
the 3-body case, as it is in the two-body case, and it may 
depend on both the family of the 3-body orbit and its 
angular momentum” [9]. Li and Liao [10] enhanced the 
aforementioned relation to T|B|>/* = T*, and numeri- 
cally proved that T* is approximately a universal con- 
stant, namely T* = 2.433 + 0.075, for the 3-body system 
with an equal mass and zero angular momentum. This 
remarkable scale-invariance period T* = T|E|?/? was a- 
gain proved by Li, Jing and Liao [11, 12] for the 3-body 
system with unequal mass and zero angular momentum, 
where T* = 3.074mz3 — 0.617 in the case of mı = mz = 1 
and mg varied. 

Although the aforementioned relations [9-11] were sup- 
ported by the numerical experiments in Ref.[9-11], The 
questions still remain whether T|E|?/? = constant is an 
universal, and if not, what form it will takes and how 
to formulate it without a further numerical simulation. 
Would we be able to find similar relations or universal 
scaling laws, for a n-body system, and, if so, how? 

The 2-body system incorporates Kepler’s three laws, 
name the law of ellipses, the law of equal areas and the 
law of harmonies. Clearly, for the 3-body system, since 
the periodic orbit is no longer elliptical, so that there is 
no corresponding law of ellipses, and because the peri- 
odic orbital topology is more complex, the law of equal 
areas might also not established. From a large number 
of numerical simulations of the 3-body system, the time 
of each object walking along its orbit is the same, that 
is, for a given mass of the 3-body system, the periodici- 
ty of the periodic orbit in the gravitational field satisfies 
Kepler’s third law, namely the law of harmonies. 

The question now is whether or not the conclusion- 
s drawn from those limited numerical experiments are 
prevalent in the 3-body system. If so, can it possibly 
be extended to a n-body system depicted in Figure 3? 
Clearly, with an increase in the number of point mass, 
the system will have more and more degrees of freedom; 
while accordingly, the dynamic process becomes hugely 
complex. If you continue to use numerical simulation to 
study the n-body problem, the calculation will certainly 
become more and more challenge. Hence, we have no 
choice but to find an alternative approach. 

This study has attempted to attack the n-body system 
(including 3-body) by using dimensional analysis [13-18]. 
The most powerful use of dimensional analysis is to pre- 
dict the outcome of an numerical experiment, depending 
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FIG. 3: n-body system 


on the variables, whilst providing theoretical insight. Di- 
mensional analysis may come across as simply trying to 
fit pieces of a puzzle together by trial and error. Howev- 
er, identifying the quantities that are relevant for a given 
problem is a demanding task, which requires deep physi- 
cal insight [18]. This may be done as follows: make a list 
of all quantities on which the answer must depend, then 
write down the dimensions of these quantities, and finally 
demand that these quantities should be combined into a 
functional form that provides the right dimension. This 
scheme was cast into a formal framework by Bucking- 
ham in 1921 and is often referred to as the Buckingham 
m-theorem [19]. 


DIMENSIONAL FORMULATION 


An important class of central gravitation fields is 
formed by those in which the potential energy is inverse- 
ly proportional to the radius r, and the force accordingly 
inversely proportional to r?. They include the fields of 
Newtonian gravitational and of Coulomb electrostatic in- 
teraction; the latter may be either attraction or repulsive. 
For the attractive gravitation field, the potential energy 
is U(r) = —a/r, where a is a positive constant with 
dimension [L]?[M][T], which is from Newtonian 2-body 
potential [3]. 

Let’s having a n-body system with N point masses, 
denoted by mz (k = 1,...,N), where each mass does 
periodic orbital motion in the Newtonian gravitational 
field, whose gravitational constant is G. Assuming that 
each point mass has no angular rotation and collision 
with each other. The question becomes how to extract 
the basic parameters of the problem from these limited 
information. The success or failure of using dimensional 
analysis depends on how to select the basic parameters 
of the problem. If one choose the wrong parameters, it 
will lead to absurd results. Hence, for the n-body system, 
which parameters are the basic parameters? 

In physics, these basic parameters should include the 
gravitational field, the characteristic area or length scale 
of the periodic orbit, as well as the characteristic mass of 
the system. It is clear that we can use the gravitational 
constant G to describe the Newtonian gravitational field; 
the characteristic mass can be the reduced mass un, Each 


cycle of the track is different, and not like the elliptical or- 
bit of the semi-long axis, which is the characteristic scale, 
but the periodic orbit has a common topological feature, 
which is that all periodic orbit are closed, hence the area 
A,, of the closed orbit could be chosen as the character- 
istic area scale of the orbit, while its square root is the 
length scale. Now the periodic orbit problem of the n- 
body system then becomes to find the orbital period, Tn, 
and total energy, |En], which is the summation of kinetic 
and potential energy, here takes its abstract value since it 
is negative for periodic orbit. From Newtonian gravita- 
tion theory, the attraction forces between bodies are lin- 
ear proportional to the product of Gm,;m,;, which means 
that the gravitation constant G can be absorbed into a 
new parameter a,,, whose dimension is [L]°[M][T]~?. In 
the following, we will use a instead of G as basic param- 
eter. 

The dimensions of those parameters are listed in the 
table below: 


TABLE I: Parameters and Dimensions 
an Hn An |Tn |En] 
LPMIT M IL? |) MIL Pe)? 


According to the dimensional analysis [13-17], the to- 
tal energy |E,,| can be expressed as 


|En| = f(Qn, Hn, An), (1) 


where f stands for a function. This relation has four 
parameters with three basic dimensions, namely time 
[T], mass [M] and length [L]. From Buckingham m- 
theorem [13-17], It can produce only one dimensionless 
parameter, t = |E,|a%u° AS, the homogenous dimen- 
sion theorem gives us a = —1, b = 0,c = 1/2, hence 
T= |En] AR Jan, which must be a constant since it is 
one only. Therefore, we have 


|En] A}? (an)! = const. (2) 


In the similar way,the orbital period Ta can be expressed 
as Tn = F(an, Mn, A), where F stands for a func- 
tion. This relation has four variables with three ba- 
sic dimensions, namely time [T], mass [M] and length 
[L]. It can produce only one dimensionless parameter, 
T= T A~a ua, which must also be a constant. 
Therefore n-body Kepler’s third law is given by 


T AR? al P uz P = const. (3) 


Combining the Eq.(2) and(3) and removing the area Ay, 
we can obtain a popular format of the Kepler’s third law 
as follows 


T,,|En|?/? = const. x an Vin. (4) 


It is noted that the period depends on the energy of the 
mass. The higher energy the system has, the shorter 


period it has. For each energy level, there is correspond- 
ing period, therefore infinite periodic orbits are existed 
(10, 11]. 


KEPLER’S THIRD LAW OF N-BODY SYSTEM 


Although Eq.(4) has been formulated, we still cannot 
get much useful information if the constant, a, and un 
cannot be defined. In other words, the success and failure 
of Eq.(4) totally depends on determination of these three 
parameters. The situation is even more difficulty since 
the only analytic information is the 2-body Kepler’s third 
law. Let’s embark our journey from Kepler’s third law. 

Since dimensional result Eq.(4) is a general result and 
should also be true for 2-body system, we can determine 
the constant by this understanding. In Eq.(4), if we can 
set the reduced mass u2 = ae and parameter ag = 
Gm,mz, by comparing with Kepler’s third law T|E|3/? = 


(mima)? ] 1"? _ 
Ta [mat] , then we can propose that const.=5. 


If we carry on and will face a challenge problem, that 
is how to extend the 2-body Kepler’s third law to n- 
body system? From Newtonian gravitation theory,for 2- 
body system (mı, M2), the attraction forces in between 
are proportional to the linear combination of m1mz; for 
3-body system (m1,m2,m3), the attraction forces be- 
tween bodies are proportional to the linear combina- 
tion of mym2, mym3 and mmg; and for 4-body sys- 
tem (m1, M2, M3, M4), the attraction forces between bod- 
ies are proportional to the linear combination of m ma, 
mM M3, MıM4, MaM3, MaM4, and m3m4; and for 5-body 
system (m1, M2, M3, M4, M5), the attraction forces be- 
tween bodies are proportional to the linear combination 
of MıM2, MIıM3Z, MıM4, MıMs5, M2M3, M2M4, M2M5, 
m3m,4, M3M5, and Mm4Ms5; and so on, for n-body, the 
attraction forces between bodies are proportional to the 
linear combination of m;m;, i = 1...n—1,j = i+1. From 
mathematics of combination, the number of combination 


of mass product is ( j= COE 


We know that Kepler’s third law 


1/2 
T2|E2|9/? = Za Mima 
V2 [m +m 


, which has only one mass product mımə of 2-body sys- 
tem. However, for 3-body system we have three mass 
product combinations, namely 


my™m2, ™ M3, M2M3 


, analogy to Kepler’s third law, let’s propose 3-body Ke- 
pler’s third law as follows 


Ty) E|?” -Tg (mimo)? + (mim3)? + (mam)? 1/2 
V2 My + mM + m3 


(5) 


Clearly, when m3 = 0, Eq.(5) is reduced to 2-body Ke- 
pler’s law. In physics, any n-body Kepler’s law should be 
able to give 2-body Kepler’s law, which means that the n- 
body Kepler’s third law must be compatible with 2-body 
Kepler’s law. In this regards, our formulation Eq.(5) is 
compatible with Kepler’s third law. 

In the light of Kepler’s third law, we would like to 
propose following conjecture: n-body Kepler’s third law 
could be expressed as follows 


where 


and total mass 


. Eq.(6) has answered the conjecture given by [6], name- 
ly, T|E|3/2 = constant, the constant on the right-hand- 
side of this equation for a specific mass system is only a 
constant rather than “universal” . 

Some examples from Eq.(6) are listed in the following 
table III. 


TABLE II: Mn and Sn of point mass system (m1, M2, ..., Mn ) 
N Mn and Sn 
2 Mz = mı + m2, S2 = (mı m2)” 
3 M3 = mı + M2 + ms, 
S3 = (mım2)? + (mims3)? + (m2m3)? 
4 Ma4 = mı + m2 + m3 + m4, 
S4= (mimz2) + (mims3 34 (mima)3 
+(mzm3)? + (mma)? + (m3ma)? 
5 Ms = mı + M2 + M3 + M4 + M5, 
Ss = (mimz2)? + (mim3)3 + (mima)? + (mıms)? 
+(m2m3)? + (mama)? + (mams)? 
+(m3ma)? + (mams)? + (mams)? 


NUMERICAL VALIDATION AND DISCUSSIONS 


To compare with numerical simulation results carried 
out by [10, 11], numerical validation will be conducted 
for following cases. 

Case 1: For 2-body system with G = 1 and mı = m2 = 
1, Kepler’s third law gives T2|E2|3/? = Z = 1.5707963. 

Case 2: For 2-body system with G = 1 and mı = 
1 and mə varied. Kepler’s third law gives T2|E2 


SUES )t/2, This law is plotted in Fig 4 for different 


range of mass mz. The law has obvious nonlinearity at 
mg < 1 and will be more linear as mz increasing. 


FIG. 4: Kepler’s law for G = 1 and mı = 1 and mə varied 


Case 3: For 3-body system with G = 1 and m = 
mz = ms = 1, Li and Liao [10] obtained T|E|°/? = 
2.433 + 0.075 through numerical curve-fitting. In this 
case, Eq.(5) produces T3|E3|3/2 = Jy = 2.22144. Our 
prediction is close to 2.358 with error 5.8%. 

Case 4: For 3-body system with G = 1 and mı = m2 = 
1 and mz varied., Li, Jing and Liao’s [11] proposed linear 
law T|E|3/? = 3.074m3 — 0.617. In this case, Eq.(5) gives 


3/25. 4 
T3|E3| 4/2 | 
This relation indicates the Kepler’s third law is nonlinear 
function of m3. The numerical comparing are illustrated 
in Fig.5 and Fig.6. 

The difference for m3 € [0,1] shown in Fig.(5), might 
be interpreted as follows: If set m3 = 0, Eq.(7) gives 
2-body Kepler’s third law T2|E2|*/? = Z; however, the 
linear law T|E|°/? = 3.074m3—0.617 obtained by Li, Jing 
and Liao [11] gives T|E|°/? = —0.617, which might have 
no physical meaning, since T|E|°/? should be positive. 
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Generally speaking, the figures show that our formulation a 
has good linearity and keep same trends as the linear law as | FE S227 seminer 
. . . . 1 ~~ 3-body system a 
T|E|3/2 = 3.074m3 — 0.617 obtained by Li, Jing and Liao fea system 
30-4 [—5-body system L 
[11]. 
If we keep mı varied and set all other point masses to Pasi [ 
be unit mass, namely mp = 1, k 41, what’s going to be 04 [ 
happen?. Four cases are listed in the below table III. 187 T 
104 H 
54 j 
TABLE III: M, and Sn of point mass system (m1, 1,1,1,1) 
N m, Mn and Sn ~~ 2 = 2 a o E 
2 my varied m2 = 1 Mz = mı +1, S2 = (m1)? m, 
3 mı varied m2 = m3 = 1 Mz = mı +2, 93 = 2(mı)? +1 
4| mı varied mp = 1, k = 2,3,4 | M4 = mı +3, S4 = 3(mi1)? +3 
5 |m; varied mp = 1, k = 2,3,4, 5| M5 = mı +4, S5 = 4(mı)” +6 FIG. 7: mı € [0,12] 


Kepler’s third law of those four cases have been illus- 
trated in the Fig.(7), which indicate that the more point 
masses the system has the higher orbit period it has. In 
general, T44|En41|?/? > Ta |En|3/ can be proven. 

For a system with point masses (M,1,1,1,..,1), if M 
were massive and much heavier than other unit mass, 


then we have interesting result as follows 
—1 
Ta |En = (—S—)'?nGM. (8) 


Clearly, it is a linear law of single massive mass M. 


CONCLUSIONS 


In summary, this study considered the periodic orbital 
period of a n-body system from the perspective of dimen- 
sional analysis. The universal law of the n-body system 
is deduced: Kepler’s third law, namely the periodic law, 
states that periodic motion of a n-body system satisfies 
the 3/4 power law of the orbital area, or the 3/2 power 
law of the total energy. In light of Kepler’s third law, 
we proposed a compatible Kepler’s third law for n-body 
system. A numerical validation and comparison study 
was hence conducted. This study may open a new venue 
for the investigation of multi-body system. 
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